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, ' ■ Abstract 

We investigate tiie stability of the Einstein static universe as a non- 
I LRS Bianchi type IX solution of the Einstein equations in the presence 

■ of both non-tilted and tilted fluids. We find that the static universe is 
^ ■ unstable to homogeneous perturbations of Bianchi type IX to the future 

^SJ and the past. 

\o- 

0\ • PACS numbers: 98.80. Cq, 98.80.Bp, 98.80. Jk 

cn: 

O ! 1 Introduction 

T— I ! 

^-H ■ The stability of the Einstein static universe is a problem with a history as old 

^ , as relativistic cosmology. The static universe was the first cosmological solution 

JJ"! ■ of general relativity to be found, by Einstein [T], and the study of its stability, 

rS I implicitly by Lemaitre [31 [3] , and then explicitly by Eddington [i] , played an 

■ important role in establishing the concept of the expanding universe. These 
authors only showed that the Einstein static solution was unstable to spatially 
homogeneous and isotropic perturbations with a limited range of perfect fluid 
equations of state. Nevertheless, after 1930, it was widely believed that the 
Einstein static universe was always unstable and would evolve into a state of 
universal contraction or expansion if perturbed [5]. In 1967, Harrison [5] showed 
that all the perturbations which render the dust-filled Einstein static universe 
unstable are in fact oscillatory if the dust is replaced by black body radiation, 
and also showed that the Einstein static universe is neutrally stable against 
small inhomogeneous perturbations for a range of perfect fluids. Similar results 
were found for homogeneous and isotropic perturbations of all amplitudes by 
Gibbons [71 [S] who also provided a thermodynamic perspective and showed that 
\i dp/dp > 1/5, where p is the pressure and p is the density of the contents of the 
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universe, then the static universe is stable. The physical reason for this surpris- 
ing conclusion is that the Einstein static space is compact and when the sound 
speed exceeds l/Vb the Jeans length exceeds the size of the universe and the in- 
stability growth is quenched by the pressure. These studies are quite limited and 
were extended by Barrow, Ellis, Maartens and Tsagas [9] to consider the whole 
spectrum of inhomogeneous scalar, vector and tensor perturbations to the static 
universe, as well as the presence of a self-interacting scalar field. They found 
that the static universe is neutrally stable against all inhomogeneous adiabatic 
scalar perturbations if dp /dp > 1/5 and confirmed that the dust {p = 0) model 
is always unstable. Any initial vector perturbations remain frozen in and so 
there is neutral stability against vortical perturbations for all equations of state 
at linear order. Likewise, there is neutral stability with respect to transverse- 
traceless tensor perturbations for all equations of state. Two caveats are to be 
noted. Since the Einstein static universe has compact space sections and many 
Killing vectors, any perturbation analysis is susceptible to the problems of lin- 
earisation instability. In effect, Einstein static is a conical point in the space 
of all solutions to the Einstein equations and special higher-order constraints 
must be checked to ensure that the leading terms in any perturbative series 
expansion around it are the leading-order terms in series that converge to true 
solutions of the Einstein equations, see for example refs [lOl [HI [9j [12]. Also, 
as an aside, we note that if 'ghost' fields with negative densities {p < 0) are 
permitted then Einstein's static universe is a stable solution of the field equa- 
tions and exact isotropic and homogeneous cosmological solutions can be found 
that exhibit stable oscillations of arbitrarily small amplitude around the static 
state |13j - in effect these are oscillating universes with non-zero minimum and 
maximum radii. 

Einstein static universes are also particular solutions in a wide range of 
extensions of general relativity. Existence conditions are known, for example, for 
gravitation theories in which the lagrangian is an arbitrary function of the scalar, 
Ricci, and Riemann curvature invariants [14j[T5]. There have been many studies 
of the Einstein static stability in other gravity theories (TT l ITS l [T9 l (20 1 [2T 1 [22 l 
[231 [23] , and from a thermodynamic point of view [3S] . There are also possible 
cosmological initial states that rely upon stabilities of the static universe, for 
example that of the emergent universe scenario [26] , which requires the Einstein 
static universe to be stable to the infinite past and expand like an Eddington- 
Lemaitre universe to the future. Eddington favoured a cosmological scenario of 
this sort, with a infinite geometrical past but only a finite thermodynamic past 
when the universe was significantly out of equilibrium. 

These earlier studies of the stability of Einstein static universes have ex- 
amined the behaviour of small and large amplitude homogeneous and isotropic 
(conformal) perturbations, and small amplitude inhomogeneous perturbations. 
However, in our earlier study of the stability properties in general relativity, 
[5], we also pointed out that a key ingredient in understanding the stability of 
Einstein static is to understand its response to long-wavelength homogeneous 
gravitational waves of Bianchi type IX. We showed that there was instability 
with respect to the mode that controlled the volume expansion in the non-tilted 
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perfect fluid case. In this paper we will provide a full analysis of the stability of 
the Einstein universe with respect to the non-LRS Bianchi type IX degrees of 
freedom for tilted and non-tilted perfect fluids. This problem reduces to study- 
ing the stability properties of the special Einstein static universe solution of 
the Bianchi type IX ('Mixmaster') Einstein equations. This anisotropic metric 
describes the most general spatially homogeneous closed universe with Einstein 
static as a special cas£0- We will show that the static universe is unstable to 
the past and to the future and we expect the solution to be unstable in other 
gravity theories as well if Mixmaster exact solutions exist. In particular, cos- 
mological scenarios with a past-eternal 'initial' Einstein static are unstable in 
general relativity. 

2 Non-tilted Bianchi Type IX Universes 

We will proceed by showing that the Einstein static universe is a particular exact 
Bianchi type IX solution of the Einstein equations for spatially homogeneous 
cosmologies and then determine its stability. This will show any effects of the 
homogeneous gravitational degrees of freedom present in type IX universes. This 
metric possesses anisotropic spatial curvature and is the only member of the 
Bianchi classification of homogeneous spaces to contain the closed Friedmann 
universes as a special case. We will examine the situation with a non-tilted 
perfect fluid source first. 

We follow Ellis and MacCallum j^7\ to write down the equations in a group 
invariant orthonormal frame. The geometry of Bianchi type IX space-time 
is characterised by the kinematic quantities of the homogeneous spatial hy- 
persurface: bulk expansion rate H, trace- free shear tensor ffap^ an auxiliary 
three-vector fla that measures the rotation of the frame with respect to Fermi- 
propagated one and a symmetric three-tensor na/s which determine the internal 
geometry of the spatial hypersurface. The spatial scalar curvature and trace- free 
part of Ricci tensor are given by 

When the fluid is non-tilted, we can diagonalise aajS and simultaneously, 
which also implies fta — 0. We assume the fluid satisfies 

P^{l-l)p (1) 

with a constant 7. The cosmological constant A can be regarded as another 
fluid with = — Pa = A. Using the notations in [55], there are eight variables 
{H, a±, ^1.2,3, P, Pa}- For Bianchi type IX, ^^l,2,3 are all non-zero and of a same 
signature, here taken to be positive. 

^Note that Einstein static is not a special case of the Kantowski-Sachs closed anisotropic 
universes, which have X S^spacc sections. 
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To obtain dimensionless equations, let us define a new variable 



D ^ <^ + i(nin2n3)3, 

which we will use to normalise the system of equations for Bianchi IX, following 
Hewitt, Uggla and Wainwright [29j. The advantage of using this variable is that 
the normalisation with respect to D is well-defined everywhere in the type IX 
state space since D 7^ as long as 711^2.3 7^ 0. Particularly, it is non-zero for 
Einstein static universe where H = 0. 

We introduce dimensionless variables as follows: 



D' D ' D' 3D^' 31)2 



We also introduce a new time variable r by defining 



dt 1 ^ 

and denote d/dr by '. 

We may use the normalised equations given in p9l with a slight modification 
to accommodate two fluids. First of all, the evolution of D decouples as 



{l + q)HD 



(3) 



where q is algebraically determined by the normalised variables through Ray- 
chaudhuri equation 
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q = 2(S: 

The rest of the evolution equations read 



S-) + ^(37-2)f^M 



H' 










-{2 - q)Hf:± - s±, 


N[ 




(i/g-4S+)7Vi, 






(i/g + 2S+ + 2\/3S_)iV2, 






{Hq + 2t+ -2V3f:_)N3, 






{2q-3j + 2)HnM, 






2{q + l)HnA. 


The components of spatial Ricci tensor are given by 




{N 


2 - Nsf - Ni{2Ni - N2 ~ m) 


2V3 L 


iN2 - NsXNi - N2 - N3) . 



(4) 



(5) 
(6) 
(7) 
(8) 

(9) 
(10) 

(11) 
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There are two constraint equations 
1 = + 

1 



v + n 



M 



where 



V 



12 



1 



(12) 
(13) 



+ 7V| + 7V| - 2N1N2 - 2N2N3 - 2N3N1 + 3 (^i^2^3) ^ 



(14) 



3 The Stability of Einstein Static Universes 

The only equilibrium point in the interior of type IX invariant set corresponds 
to the Einstein static solution, given by 



F = E+ = 0, 



^1,2.3 = 2, 



37 



r^A = 1 ■ 

37 



(15) 



We are interested in linear stability around it. 
Linearisation: 

We denote a small deviation of a quantity x from its value at the Einstein 
static solution, xq, hy Sx — x ^ xq. From (jl3p . we have 

(16) 
(17) 



SNi + SN2 + SN3 = 0, 
and, with ([T^ . this means that 

SV = 0, SQm = -StlA. 
From the definitions of q and 5*^, we have 

1 



Sq 

SS+ 
5S^ 



-l{2SNi - 6N2 -SN^)^ SN2 + SN3 
3 



V3 



{SN3 -dN2). 



(18) 
(19) 
(20) 



We replace SNi,2,3 by SS± as the variables for the linearised equations. Note 
that the relation (IT51) was used to reduce the number of variables. Finally, the 
evolution of the perturbations are given by 



6H' 




SQ'a 


..(.- 




= -6S± 


SS'± 


= 8(5E±. 



(21) 
(22) 

(23) 
(24) 
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The linear stability of the static solution is therefore decided by the eigenvalues: 

A = ±V37- 2, ±i2V2, ±i2y/2. (25) 

We see that there is always an instability in the direction of increasing and de- 
creasing T time for 7 > 2/3 because of the two real eigenvalues with opposite 
signs. Note that this mode is associated with the volume expansion (it is the 
one identified in ref [9]), namely the FLRW mode and would still be present 
in the isotropic limit where iVi 2,3 — Sj- — 0. The effects of the anisotropic 
curvature are not decisive at this order and produce the two pairs of purely 
imaginary eigenvalues. We note that two of the imaginary eigenvalues corre- 
spond to an axisymmetric mode (5S+ — ^S^) mentioned by ^D] who studied 
only the locally rotationally symmetric (LRS) case, although their analysis con- 
sidered occurrence of chaos [31] following evolution away from the Einstein static 
universe. 

4 The Extension to Tilted Fluids 

We can add a tilt to the fluid velocity [33], so that the unit normal of the 
homogeneous hypersurfaces and the fluid 4-velocity u°' are not coincident, 
but differ by the addition of a peculiar velocity vector v°- , with associated scalar 

11, so that 

^vin" ^ v") r = (i-w2)-i 

Denote the energy density of the fluid seen by the observer by p and it satisfies 
the equation of state H]). If we decompose its energy-momentum tensor with 
respect to as 

Tah = fJ-nanb + 2q(anb) + P{9ab + "a^b) + T^ab 

such that qaTi"" = tt^^ — and iTab'n'' = 0, then these quantities are given by 

^l = (i + 7rV)p, (26) 

p - _ 1 + p, (27) 

Qa = iT^pVa, (28) 



71"ab = 7r 



VaVb - ^v'^igab + naUb) 



(29) 



In this case, the orthogonal equations in [29] are no longer usable and we start 
from the general equations for an arbitrary group invariant orthonormal frame 
given in [32j . Since the fluid velocity gives rise to an energy flux with respect 
to the observer n", from the equations (1.93) in [32], we cannot diagonalise a a/3 
and nai3 simultaneously. Although we could use the remaining gauge freedom 
to choose the spatial triad to diagonalise Uap, this does not turn out to be the 
best gauge choice for analyzing the stability of isotopic equilibrium points, as 
we will see. Instead, we proceed without fixing the gauge, which means that we 
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need to modify the definition of the normalisation factor D. Noting that for an 
arbitrary matrix A, 

^detA = detAtr(^A-i^^ , (30) 

it is easy to see from the equation (1.96) in [32 that 

^(detn„/3) = detn,;3 

= — 3i7detnQ^ 

and so if we redefine D as 



D= -^/i/2 + i(detn„0)^ 



then the normaUsation wih work as before. Let us first define the normafised 
geometrical quantities by 

There now exists an angular velocity, fJo,, of the frame with respect to the 
Fermi-propagated frame. This is non-dynamical and is just a manifestation of 
the gauge freedom. It is normalised by 

Ka = 

D 

The fluid density parameter seen by the observer n° is defined by 

Since the tilt velocity, is already dimensionless, it will not be normalised. 
The cosmological constant represents another dynamical degree of freedom through 
the variable 

3D^ 

We can now derive a new system of equations from the equations (1.90)-(1.100) 
in [32]. 

First of all, the evolution of D decouples as before; 

b' = -{l + q)H^b (31) 
where ' now denotes d/dr ~ D^^d/dt and q is determined by 



q = 2f:-i:,. %l\^l-^^:^ ^M - h,. (32) 
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The rest of the Einstein equations read 



1 = ±^''±^,+V + nM+^A, 







1 _(^_ 1)^,2 %^V/3,, 



where enua is the spatial Levi-Civita symbol and 



Nap - ,:N^^5ap 



V 



3717 



M 



1 - (7-1)^2 

1 / . \ i 1 

4V "^y '"6 

The Jacobi identities become 



aP 



N'ap = qHNap + 2t^^Np)^ + 2e^'\^Np)^R,. 

The evolution for the matter variables is obtained by substituting ((26 
the contracted Bianchi identities: 



(33) 
(34) 

(35) 



(36) 
(37) 
(38) 

(39) 
29t into 



M 



2{l + q)HnM- 
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1 + (7- l)i;2 

r-2(37 - 4)g - (7 - 2)±^,v^v'' 
1 - (7 - 1)^2 

-Sq^w'^ - eaf_ii,R''v'' + eativN^i^v'^v" , 

2{l + q)HnA 



{3 + v^)H + ±apv"v^ 



Finally, by the definition of the normalisation, we have 

H' = -{l-H^)q, 

1 ~ 

1 = H^ + -{detNapy\ 



^M, (40) 
(41) 

(42) 

(43) 
(44) 



The normalised fluid density, ^Im, is usually eliminated by the constraint 
(|34|) and then there are 16 dynamical variables (taking into account the traceless 
condition for Eq^). We still have 3 degrees of gauge freedom and 4 constraints 
( (|35|) and (jS])), so the dynamical system is 9-dimensional, which is consistent 
with the counting of the ref [34j because now we have one more degree of freedom 
coming from addition of the cosmological constant to the 8 in general single- fluid 
tilted models. 

To find equilibrium points in the system, we need consider the gauge. Unless 
we exhaust this freedom, an equilibrium point can be time dependent. A typical 
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choice is to diagonalize Nafi- It leads to the following defining equations for the 
non-dynamical Re,: 

Ri 
R2 
R3 



S23(-^22 


+ N33) 


N22 - 


N33 


^13 (-^33 


+ ^11) 


N33- 


Nn 


Sl2(^ll 


+ N22) 


iVn- 


N22 



If we substitute those expressions into the evolution equations, they become 
singular when the denominators vanish, which includes the isotropic case. Thus 
this gauge is not suited for the stability analysis of Einstein static universe. The 
situation is the same for a gauge where E^^ is diagonal. 

Fortunately, in the Einstein static universe all the dynamical variables have 
to be time independent since all the gauge dependent variables are zero. While 
there is ambiguity in the values of Ra which are arbitrary (potentially time- 
dependent) since they merely represent the rotation of the spatial frame aside 
from it, the equilibrium point is gauge invariantly characterised by 

2 

Let us consider a perturbation around this background without fixing the gauge. 
In this way, we expect to have 12 eigenvalues taking into account the four 
constraints mentioned above. From now on, the quantities without S prefixes 
are understood to take the background values. The constraint ((35)) becomes 

Thus, the velocity perturbation to linear order automatically vanishes, and so 
6Ilap — 0. Secondly, the constraint (j44|) implies 

= (5(detiV„/j) 
d det Nai3 



where we used ([50]) again. Noting that 



and 



dN 
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we can see that 

<5iV/ - 0, 

namely, SNap is trace-free. This ensures that SV = 0. It's also easy to see that 
SSap = 26 Nap. Therefore, the linearised evolution equations are quite simple: 

S±'ap - ^2SNap+2e''\jtp)^K, 
SN'^p = AStap + 2e'"'^JNp)^R,, 

SCl'^ = 2(^l-^yH. 

The equations decouple into two. The scalar part is exactly the same as the 
non-tilted case and gives one positive and one negative eigenvalue and so the 
Einstein static universe is also unstable under the same conditions as in the non- 
tilted case. Therefore there is instability of the static solution for increasing and 
decreasing times. The tensorial part can be brought into a standard form 

St'ap = -25Nafi 

by an appropriate rotation of the spatial frame, leaving five pairs of coupled 
equations, and therefore five sets of eigenvalues ±i2y/2. Although not all of 
them are physical as we haven't fix the gauge, this is unimportant because all 
the eigenvalues are purely imaginary. This means the nature of the tensorial 
perturbations is also the same as in the non-tilted case. 



5 Conclusions 

In this study of the stability of the Einstein static universe against Bianchi type 
IX modes we have extended the earlier studies of ref [9] , that considered only the 
destabilizing efl[ects of a single homogeneous mode in the presence of comoving 
fiuids, and those of ref [W, who considered only the evolution of LRS type IX 
models with non-tilted fluid. We investigated the general type IX evolution in 
the vicinity of the Einstein static model in the presence of a fluid with non- 
tilted and tilted fluid motion. We have established that the Einstein static 
universe is unstable to Bianchi type-IX spatially homogeneous perturbations 
in the presence of non-tilted and tilted perfect fluids with p + 3p > 0. We 
also found that the imaginary eigenvalues corresponding to the perturbative 
effects of anisotropic curvature (Mixmaster modes) and fluid tilt generalise the 
oscillatory behaviour of the finite wavelength vector and tensor perturbations 
found in early studies of small amplitude perturbations. 
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